Alternative RNA splicing defects in pediatric cancers: new insights in tumorigenesis and potential therapeutic vulnerabilities.

BACKGROUND: Compared with adult cancers, pediatric cancers are uniquely characterized by a genomically stable landscape and lower tumor mutational burden. Alternative splicing, however, a global cellular process that produces different messenger RNA/protein isoforms from a single messenger RNA transcript, has been increasingly implicated in the development of pediatric cancers. DESIGN: We review the current literature on the role of alternative splicing in adult cancer, cancer predisposition syndromes, and pediatric cancers. We also describe multiple splice variants identified in adult cancers and confirmed through comprehensive genomic profiling in our institutional cohort of rare, refractory, and relapsed pediatric and adolescent young adult cancer patients. Finally, we summarize the contributions of alternative splicing events to neoantigens and chemoresistance and prospects for splicing-based therapies. RESULTS: Published dysregulated splicing events can be categorized as exon inclusion, exon exclusion, splicing factor up-regulation, or splice site alterations. We observe these phenomena in cancer predisposition syndromes (Lynch syndrome, Li-Fraumeni syndrome, CHEK2) and pediatric leukemia (B-cell acute lymphoblastic leukemia), sarcomas (Ewing sarcoma, rhabdomyosarcoma, osteosarcoma), retinoblastoma, Wilms' tumor, and neuroblastoma. Within our institutional cohort, we demonstrate splice variants in key regulatory genes (CHEK2, TP53, PIK3R1, MDM2, KDM6A, NF1) that resulted in exon exclusion or splice site alterations, which were predicted to impact functional protein expression and promote tumorigenesis. Differentially spliced isoforms and splicing proteins also impact neoantigen creation and treatment resistance, such as imatinib or glucocorticoid regimens. Additionally, splice-altering strategies with the potential to change the therapeutic landscape of pediatric cancers include antisense oligonucleotides, adeno-associated virus gene transfers, and small molecule inhibitors. CONCLUSIONS: Alternative splicing plays a critical role in the formation and growth of pediatric cancers, and our institutional cohort confirms and highlights the broad spectrum of affected genes in a variety of cancers. Further studies that elucidate the mechanisms of disease-inducing splicing events will contribute toward the development of novel therapeutics.

Genomic characterization of HER2-positive breast cancer and response to neoadjuvant trastuzumab and chemotherapy-results from the ACOSOG Z1041 (Alliance) trial.

Background: HER2 (ERBB2) gene amplification and its corresponding overexpression are present in 15-30% of invasive breast cancers. While HER2-targeted agents are effective treatments, resistance remains a major cause of death. The American College of Surgeons Oncology Group Z1041 trial (NCT00513292) was designed to compare the pathologic complete response (pCR) rate of distinct regimens of neoadjuvant chemotherapy and trastuzumab, but ultimately identified no difference. Patients and methods: In supplement to tissues from 37 Z1041 cases, 11 similarly treated cases were obtained from a single institution study (NCT00353483). We have extracted genomic DNA from both pre-treatment tumor biopsies and blood of these 48 cases, and performed whole genome (WGS) and exome sequencing. Coincident with these efforts, we have generated RNA-seq profiles from 42 of the tumor biopsies. Among patients in this cohort, 24 (50%) achieved a pCR. Results: We have characterized the genomic landscape of HER2-positive breast cancer and investigated associations between genomic features and pCR. Cases assigned to the HER2-enriched subtype by RNA-seq analysis were more likely to achieve a pCR compared to the luminal, basal-like, or normal-like subtypes (19/27 versus 3/15; P = 0.0032). Mutational events led to the generation of putatively active neoantigens, but were overall not associated with pCR. ERBB2 and GRB7 were the genes most commonly observed in fusion events, and genomic copy number analysis of the ERBB2 locus indicated that cases with either no observable or low-level ERBB2 amplification were less likely to achieve a pCR (7/8 versus 17/40; P = 0.048). Moreover, among cases that achieved a pCR, tumors consistently expressed immune signatures that may contribute to therapeutic response. Conclusion: The identification of these features suggests that it may be possible to predict, at the time of diagnosis, those HER2-positive breast cancer patients who will not respond to treatment with chemotherapy and trastuzumab. ClinicalTrials.gov identifiers: NCT00513292, NCT00353483.

CTCF genetic alterations in endometrial carcinoma are pro-tumorigenic.

CTCF is a haploinsufficient tumour suppressor gene with diverse normal functions in genome structure and gene regulation. However the mechanism by which CTCF haploinsufficiency contributes to cancer development is not well understood. CTCF is frequently mutated in endometrial cancer. Here we show that most CTCF mutations effectively result in CTCF haploinsufficiency through nonsense-mediated decay of mutant transcripts, or loss-of-function missense mutation. Conversely, we identified a recurrent CTCF mutation K365T, which alters a DNA binding residue, and acts as a gain-of-function mutation enhancing cell survival. CTCF genetic deletion occurs predominantly in poor prognosis serous subtype tumours, and this genetic deletion is associated with poor overall survival. In addition, we have shown that CTCF haploinsufficiency also occurs in poor prognosis endometrial clear cell carcinomas and has some association with endometrial cancer relapse and metastasis. Using shRNA targeting CTCF to recapitulate CTCF haploinsufficiency, we have identified a novel role for CTCF in the regulation of cellular polarity of endometrial glandular epithelium. Overall, we have identified two novel pro-tumorigenic roles (promoting cell survival and altering cell polarity) for genetic alterations of CTCF in endometrial cancer.

Consensus on precision medicine for metastatic cancers: a report from the MAP conference.

Recent advances in biotechnologies have led to the development of multiplex genomic and proteomic analyses for clinical use. Nevertheless, guidelines are currently lacking to determine which molecular assays should be implemented in metastatic cancers. The first MAP conference was dedicated to exploring the use of genomics to better select therapies in the treatment of metastatic cancers. Sixteen consensus items were covered. There was a consensus that new technologies like next-generation sequencing of tumors and ddPCR on circulating free DNA have convincing analytical validity. Further work needs to be undertaken to establish the clinical utility of liquid biopsies and the added clinical value of expanding from individual gene tests into large gene panels. Experts agreed that standardized bioinformatics methods for biological interpretation of genomic data are needed and that precision medicine trials should be stratified based on the level of evidence available for the genomic alterations identified.

A genomic case study of mixed fibrolamellar hepatocellular carcinoma.

BACKGROUND: Mixed fibrolamellar hepatocellular carcinoma (mFL-HCC) is a rare liver tumor defined by the presence of both pure FL-HCC and conventional HCC components, represents up to 25% of cases of FL-HCC, and has been associated with worse prognosis. Recent genomic characterization of pure FL-HCC identified a highly recurrent transcript fusion (DNAJB1:PRKACA) not found in conventional HCC. PATIENTS AND METHODS: We performed exome and transcriptome sequencing of a case of mFL-HCC. A novel BAC-capture approach was developed to identify a 400 kb deletion as the underlying genomic mechanism for a DNAJB1:PRKACA fusion in this case. A sensitive Nanostring Elements assay was used to screen for this transcript fusion in a second case of mFL-HCC, 112 additional HCC samples and 44 adjacent non-tumor liver samples. RESULTS: We report the first comprehensive genomic analysis of a case of mFL-HCC. No common HCC-associated mutations were identified. The very low mutation rate of this case, large number of mostly single-copy, long-range copy number variants, and high expression of ERBB2 were more consistent with previous reports of pure FL-HCC than conventional HCC. In particular, the DNAJB1:PRKACA fusion transcript specifically associated with pure FL-HCC was detected at very high expression levels. Subsequent analysis revealed the presence of this fusion in all primary and metastatic samples, including those with mixed or conventional HCC pathology. A second case of mFL-HCC confirmed our finding that the fusion was detectable in conventional components. An expanded screen identified a third case of fusion-positive HCC, which upon review, also had both conventional and fibrolamellar features. This screen confirmed the absence of the fusion in all conventional HCC and adjacent non-tumor liver samples. CONCLUSION: These results indicate that mFL-HCC is similar to pure FL-HCC at the genomic level and the DNAJB1:PRKACA fusion can be used as a diagnostic tool for both pure and mFL-HCC.

Recurrent DNMT3A mutations in patients with myelodysplastic syndromes.

Alterations in DNA methylation have been implicated in the pathogenesis of myelodysplastic syndromes (MDS), although the underlying mechanism remains largely unknown. Methylation of CpG dinucleotides is mediated by DNA methyltransferases, including DNMT1, DNMT3A and DNMT3B. DNMT3A mutations have recently been reported in patients with de novo acute myeloid leukemia (AML), providing a rationale for examining the status of DNMT3A in MDS samples. In this study, we report the frequency of DNMT3A mutations in patients with de novo MDS, and their association with secondary AML. We sequenced all coding exons of DNMT3A using DNA from bone marrow and paired normal cells from 150 patients with MDS and identified 13 heterozygous mutations with predicted translational consequences in 12/150 patients (8.0%). Amino acid R882, located in the methyltransferase domain of DNMT3A, was the most common mutation site, accounting for 4/13 mutations. DNMT3A mutations were expressed in the majority of cells in all tested mutant samples regardless of myeloblast counts, suggesting that DNMT3A mutations occur early in the course of MDS. Patients with DNMT3A mutations had worse overall survival compared with patients without DNMT3A mutations (P=0.005) and more rapid progression to AML (P=0.007), suggesting that DNMT3A mutation status may have prognostic value in de novo MDS.

A map of human genome sequence variation containing 1.42 million single nucleotide polymorphisms.

We describe a map of 1.42 million single nucleotide polymorphisms (SNPs) distributed throughout the human genome, providing an average density on available sequence of one SNP every 1.9 kilobases. These SNPs were primarily discovered by two projects: The SNP Consortium and the analysis of clone overlaps by the International Human Genome Sequencing Consortium. The map integrates all publicly available SNPs with described genes and other genomic features. We estimate that 60,000 SNPs fall within exon (coding and untranslated regions), and 85% of exons are within 5 kb of the nearest SNP. Nucleotide diversity varies greatly across the genome, in a manner broadly consistent with a standard population genetic model of human history. This high-density SNP map provides a public resource for defining haplotype variation across the genome, and should help to identify biomedically important genes for diagnosis and therapy.

A physical map of the human genome.

The human genome is by far the largest genome to be sequenced, and its size and complexity present many challenges for sequence assembly. The International Human Genome Sequencing Consortium constructed a map of the whole genome to enable the selection of clones for sequencing and for the accurate assembly of the genome sequence. Here we report the construction of the whole-genome bacterial artificial chromosome (BAC) map and its integration with previous landmark maps and information from mapping efforts focused on specific chromosomal regions. We also describe the integration of sequence data with the map.

Association between DQB1 and cervical cancer in patients with human papillomavirus and family controls.

OBJECTIVE: The role of human leukocyte antigen (HLA) DQB1 alleles and human papillomavirus (HPV) as contributing factors to invasive cervical cancer was investigated. To overcome problems of misleading causal inferences common in traditional case-control studies, a family-based test, the transmission/disequilibrium test, was used. METHODS: Ninety-six patients with pathologically confirmed invasive cervical cancer were ascertained. Human papillomavirus types were determined in 80 patients, of whom 81.25% were HPV-positive, and 18.75% were HPV-negative. Deoxyribonucleic acid was extracted from samples, taken from patients and their parents, and sequenced to determine DQB1 genotypes. Nuclear family data were used to test whether the DQB1 locus is associated with invasive cervical cancer while controlling for high-risk HPV-positive patients. The transmission/disequilibrium test evaluates whether the frequency of transmission of parental marker alleles to their affected offspring deviates from the expected Mendelian frequency of 50%. RESULTS: The HLA DQB1 locus showed evidence for allelic association with invasive cervical cancer in high-risk HPV-positive patients (P = .006). The transmission/disequilibrium test showed that the DQB1*0303 allele was transmitted to high-risk HPV patients more often than expected by chance, chi2(1) = 8.0, P = .005 (P = .035 when correcting for multiple tests). Tests of association were negative when applied to all 96 patients, irrespective of HPV status. No significant differences were found in the distribution of the DQB1 alleles among HPV-positive patients compared with those who were HPV-negative, indicating that HLA alleles are not associated with susceptibility to HPV infection. CONCLUSION: These results suggest that the DQB1*0303 allele increases the risk for invasive cervical cancer in women who are HPV-positive.

Capillary electrophoresis platforms for DNA sequence analysis.

Commercial availability of multicapillary DNA sequencing instruments can significantly impact the existing paradigm in high-throughput DNA sequencing facilities, shifting the rate-limiting step away from the electrophoresis and detection step. These instruments, although novel and poorly understood, represent a significant capital investment.This review informs the reader about the state-of-the-art instruments, although capillary sequencers are in a state of almost constant flux, with improvements emerging rapidly.

A pneumatic device for rapid loading of DNA sequencing gels.

This work describes the design and construction of a device that facilitates the loading of DNA samples onto polyacrylamide gels for detection in the Perkin Elmer/Applied Biosystems (PE/ABI) 373 and 377 DNA sequencing instruments. The device is mounted onto the existing gel cassettes and makes the process of loading high-density gels less cumbersome while the associated time and errors are reduced. The principle of operation includes the simultaneous transfer of the entire batch of samples, in which a spring-loaded air cylinder generates positive pressure and flexible silica capillaries transfer the samples. A retractable capillary array carrier allows the delivery ends of the capillaries to be held up clear of the gel during loader attachment on the gel plates, while enabling their insertion in the gel wells once the device is securely mounted. Gel-loading devices capable of simultaneously transferring 72 samples onto the PE/ABI 373 and 377 are currently being used in our production sequencing groups while a 96-sample transfer prototype undergoes testing.

End sequence determination from large insert clones using energy transfer fluorescent primers.

Genome mapping strategies depend heavily on confirmatory data of several types to establish overlaps between contiguous stretches of cloned DNA derived from genomic regions. One type of ancillary data that can contribute to establishing these overlaps is DNA sequence data derived from the ends of large (> 30 kb) inserts in genomic clones. This type of data can be difficult to obtain routinely, because large clones are often unstable and microgram quantities of highly purified DNA are required in each sequencing reaction to obtain sufficient signal for accurate base calling and maximum read length. Recently, we have been experimenting with methods to consistently obtain up to 800 bases of high-quality sequence data from the ends of large insert clones using ThermoSequenase DNA polymerase and Energy Transfer fluorescent primers. Our experimental approach and results, described in this paper, indicate that routinely obtaining high-quality sequence data from the ends of large insert genomic clones is feasible. Such data can contribute to the assessment of common regions between large insert clones, to the establishment of conservation of synteny between closely related species, and to the detection of additional contiguous clones.

An infrared fluorescent dATP for labeling DNA.

Near-infrared fluorescence provides a nonradioactive method of detection with high sensitivity and low background. An infrared fluorophore has been attached covalently to the nucleotide deoxyadenosine triphosphate (dATP) to provide a reagent for enzymatic labeling of various types of DNA molecules and for facilitating their detection with an automated DNA sequencing and analysis system. DNA sequencing reaction products can be labeled internally by performing limited polymerization utilizing infrared-labeled dATP (IR-dATP) as the sole source of adenine deoxynucleotide prior to a dideoxy-specific termination reaction. PCR products can be labeled fluorescently by the addition of limited quantities of IR-dATP to the amplification reaction. This latter strategy has been utilized for detection of short tandem repeat polymorphisms (STRPs) which are useful for gene mapping, genetic diagnostics, forensic analysis, and paternity testing. Restriction fragments can be labeled also by fill-in reactions of appropriate 5' overhangs. Diminutive amounts of such fluorescently labeled DNA molecules can be visualized rapidly and conveniently using infrared detection technology.

Prep-A-Gene: a superior matrix for the purification of DNA and DNA fragments.

A new chromatographic matrix, Prep-A-Gene, is described for the isolation and purification of high quality DNA suitable for restriction analysis, ligation, transformation and sequencing protocols. This matrix selectively binds DNA greater than approximately 200 base pairs in length, while RNA, proteins, cellular components, agarose and other contaminants are washed free in minutes. This eliminates the need for time-consuming and laborious RNase treatments, gel extractions and phenol extractions. The DNA that is desorbed from the matrix is available immediately as a substrate for subsequent protocols. DNA purified in this manner exhibits no detectable shearing, even with more fragile chromosomal DNA.

Automated methods for single-stranded DNA isolation and dideoxynucleotide DNA sequencing reactions on a robotic workstation.

Automated procedures have been developed for both the simultaneous isolation of 96 single-stranded M13 chimeric template DNAs in less than two hours, and for simultaneously pipetting 24 dideoxynucleotide sequencing reactions on a commercially available laboratory workstation. The DNA sequencing results obtained by either radiolabeled or fluorescent methods are consistent with the premise that automation of these portions of DNA sequencing projects will improve the reproducibility of the DNA isolation and the procedures for these normally labor-intensive steps provides an approach for rapid acquisition of large amounts of high quality, reproducible DNA sequence data.